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1. Introduction 

Inspiralling compact binary systems are among the most promising sources for the 
emission of gravitational waves detectable with present day's gravitational wave 
interferometers. Earth bound gravitational wave detectors such as Geo600, VIRGO, 
TAMA and LIGO are most sensitive at wavelengths of about 10-1000 Hz. This 
corresponds roughly to the last 10 minutes of the inspiral before final plounge. In 
that regime Newtonian mechanics is not valid and post-Newtonian approximation must 
be applied. An earlier stage of the inspiral process will be covered by the yet to be 
built LISA interferometer. It is expected that LISA will be sensible to gravitational 
waves with frequencies from about 10~^ to 10""^ Hz. However, in order to actually 
detect gravitational waves highly accurate templates are essentially. Post-Newtonian 
corrections must be included into the EoM. Moreover, computing the gravitational 
waveforms beyond the leading order approximation higher mass and current multipole 
moments must be taken into account. 

In the past the study of close compact binaries was often based on the assumption that 
the stars can be treated as pointlike, non-spinning objects. Upon this assumption it is 
possible to derive an analytic, so called quasi-Keplerian solution for the conservative part 
of the EoM up to the third post-Newtonian approximation [1]. Dissipative effects due to 
the emission of gravitational waves first appear at the order (f/c)^, which corresponds 
to the 2.5 post-Newtonian order. 

In the presence of at least one spinning component a post-Newtonian spin-orbit coupling 
which first appears at 1.5 post-Newtonian order leads to complifications, which hamper 
the investigation of spinning compact binary systems enormeously. In general, neither 
orbital angular momentum nor the stellar spin are conserved. Up to now an analytic 
solution to the conservative part of the EoM including spin-orbit coupling has been 
found in a few special cases only (see e.g. [2]). 

Thus far not much progress has been made in the investigation of close binary systems of 
finite size objects. Within the framework of post-Newtonian analysis it is often argued 
that finite size effects are neglegible during most of the inspiral process and will become 
important not until the last few orbits before the final plunge [3]. Growing interest in 
the role of finite size effects comes mainly from the side of numerical relativity. In fact, 
though finite size effects due to stellar rotation and oscillation in compact binary systems 
are very small they can well be in the order of the first post-Newtonian corrections to 
the orbital dynamics. These secular effects, being of Newtonian origin, accumulate over 
a large number of orbits and thus, seen at longer terms, lead to significant phase shifts in 
the gravitational waves. The influence of stellar oscillations on the dynamical evolution 
and leading order gravitational wave emission has been investigated by Kokkotas and 
Schafer and Lai and Ho [5] for nonrotating, polytropic neutron stars and by Lai et 
al [6] and Hansen [7] for Riemann-S binaries. In all these approaches the analysis was 
based on Newtonian theory, the 2.5pN radiation reaction terms being the only post- 
Newtonian terms included. 
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This paper is devoted to the investigation of the influence of finite size effects on the 
dynamics and gravitational wave emission beyond the leading order approximation. 
Basically, perturbations of the point particle dynamics arise due to stellar oscillations 
(mainly tidally driven) and rotational deformation. Here we shall restrict ourselves to 
systems were the apsidal motion due to rotational deformation is much larger than the 
one caused by stellar oscillations. 

To begin with, let us note that there is, strictly speaking, no spinning point particle. A 
spinning object automatically gains a finite size (see e.g. p]). It hat been long known 
that the coupling of the nonvanishing stellar quadrupole moment to the orbital motion, 
which is a purely Newtonian effect, leads to an apsidal motion. For a couple of close 
main sequence star binaries this apsidal motion has been observed to great accuracy 
(see e.g. Claret and Willems [9]). In these systems apsidal motion due to finite size 
effects is considerably larger than the relativistic periastron advance. Another group 
of binary systems, where finite size effects play an important role, are binary pulsars 
such as PSR B 1259-63, which has been found by Johnston et al. in 1992 [10]. While 
the compact component of PSR 1259-63 is a 47 ms pulsar, it's companion is a Be star, 
whose spin-induced quadrupole deformation leads to an apsidal motion. 
However, it is close compact binary systems which are most relevant for gravitational 
wave detectors. This includes not only black hole-black hole, black hole-neutron star 
or neutron star- neutron star binaries but, if LISA is ready to work, also close white 
dwarf binaries. Of course the rotational deformation of compact stars is much smaller 
than it would be possible in non-compact stars. However, at least for neutron star 
(NS) and white dwarf (WD) binaries one should not a priori neglect finite size effects 
due to quadrupole deformation. In particular it is well possible that the perturbations 
introduced by the coupling of the stellar quadrupole moment to the orbital motion is of 
the same order of magnitude as the first post-Newtonian corrections. This is assumed 
throughout the paper. The analysis applies to compact binary systems, whose spinning 
component is a WD or a fast rotating NS. In order to simplify calculations it is assumed 
that the spin is perpendicular to the orbital plane. In section [2] the orbital evolution of a 
spinning compact binary is studied up to first post Newtonian order in the point particle 
dynamics. The EoM as well as a parametric, quasi-Keplerian solution are derived. In 
section [3] the next to leading order gravitational waveforms are calculated explicitly. 
The long time evolution and the influence of the quadrupole coupling to the inspiral 
process is discussed in section ID 

2. The IpN orbital motion including spin effects due to rotational 
deformation 

In 1985 Damour and Deruelle [11] succeeded in deriving an analytic solution to the 
IpN EoM of a point-particle binary, which exhibits a remarkable similarity to the well 
known Kepler parametrization in Newtonian theory and is hence called quasi-Keplerian 
solution. Using the same strategy Wex considered a binary consisting of a pulsar and 
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a spinning main sequence component at Newtonian order [12]. Treating the Newtonian 
couphng between the rotationally deformed star and the orbital dynamics as a small 
perturbation he derived a quasi-Keplerian solution up to first order in the deformation 
parameter g, which will be introduced in the following. In this letter we shall extend 
his investigations, taking into account the orbital dynamics up to first post-Newtonian 
approximation. In order to derive an analytic solution we shall further assume that 
the modifications induced by finite size effects are of the same order as the IpN orbital 
corrections. That is, we restrict our analysis to compact binary systems consisting of 
a fast spinning neutron star or white dwarf and a non-spinning compact object. Of 
course, all results can be applied to a mean sequence star-compact star binary in the 
Newtonian limit. 

The rotational deformation of a spinning star of mass m can be described by some 
parameter g, which is defined as [13j 

mg := ^ j dV'p{r') [r'^ - 3{s ■ r'f] = A/, (1) 

where m = J p{r')dV' is the stellar mass and A J is the difference of the moments of 
inertia parallel and perpendicular to the spin axis s, respectively. In particular, for 
rotating fluids, q is given by (see e.g. [H]) 

q = -kR^d\ Q = ^ (2) 
3 ^Gm/R^ ^ ' 

Here R denotes the polar radius and Vt the angular velocity of the rotating star. The 
constant of apsidal motion k strongly depends on the density distribution. It vanishes 
if all mass is concentrated in the center and takes it's maximal value kmax = 0.75 for a 
homogeneous sphere. Neutron stars and low mass white dwarfs can be approximately 
modelled by an polytropic EoS with index n=0.5 ... 1 and n ^ 1.5 for neutron stars and 
white dwarfs, respectively. In fact, for polytropes the value of k is given by (see e.g. 

m) 

k = l(A2(n) - 1), 

where the function A2{n) has been introduced by Chandrasekhar [T6] . 
The quadrupole deformation gives rise to a quadrupole coupling, which modifies the 
orbital dynamics as well as the gravitational wave emission of the binary. Neglecting 
contributions arising from tidally induced stellar oscillations the quadrupole coupling 
can be described by a Hamiltonian 

GMpq 



^1 2r3 



^ (3) 



where r is the orbital separation, while and fi denote total and reduced mass, 
respectively. As we shall see in the following this contribution will lead to a periastron 
advance already at Newtonian order. 

In general, the direction of the star's spin won't be conserved, which complicates the 
analysis enormeously. However, one might assume scenarios where the spin is parallel 
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to the angular momentum of the system, i.e. perpendicular to the orbital plane. In 
these cases Eq. ([3]) takes a rather simple form, 

- (4) 

and we shall restrict to this special situation further on. 

Before we proceed it is important to compare the perturbation introduced by quadrupole 
coupling with the first post-Newtonian correction terms to the EoM. Let us assume a 
spinning binary system in a circular orbit. Not taking into account the g-coupling the 
orbital energy read^ 

_ GM^^ i-vq^M^_ 

" + ~S r^ = + 

Comparing £^pjv with the coupling energy we find that the g-term offers a contribution 
comparable to the IpN orbital perturbation if 
q GM 



rs- 



If q/rQ is much larger than the Schwarzschild radius rs the Newtonian quadrupole 
contribution will clearly dominate, while for q/vo -C GM./(? the leading perturbation 
to the Keplerian orbit comes from the post-Newtonian correction terms. The value of q 
crucially depends on the density distribution (via k) and on the angular velorcity of the 
rotating star. The later one is bounded by the mass-shedding limit. For a Newtonian 
star with a polytropic EoS one can show that the mass-shedding limit is given by [17] 



where R is the polar radius of the star. Thus, for polytropic stars q is bounded by 

For a typical neutron star with polytropic index n = 1 and polar radius i? = 10 km 
the maximal value of q is qmax = 5.1 km^, while for a white dwarf with n = 1.5 
and a radius of 1000 km this value is given bj|§| qmax = 28563 km^. Now the post- 
Newtonian approximation is valid only outside the innermost stable circular orbitljj]. 
For finite size binaries the orbital separation should be considerably larger than this 
value. To be more precisely, since we do not consider mass overflow the spinning star 
should remain well inside it's Roche volume throughout this calculation. Taking this 
into account it becomes clear that for close NS-BH and NS-NS binaries the contribution 
of the quadrupole coupling is by a factor 100 or more smaller than the first post- 
Newtonian correction. However, at least for sufficiently fast rotation the g-term gives 

t Note, that £ = H. 

§ For n — 1 the constant of apsidal motion is fc = 0.259925, for n = 1.5 one finds k = 0.1446. 
II In isotrope coordinates the radius of the last stable orbit of a test particle orbiting a Schwarzschild BH 
is given by rjsco = 5GA4/c^. For compact binaries a calculation of rjsco in full General Relativity 
has not been successfuU until now, but post-Newtonian analysis suggests a value around 5GA4/c^. 
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rise to perturbations which are considerably larger than the 1.5 pN order. For spinning 
white dwarfs things can be different. In that case the finite size contribution due to 
rotational deformation can be equal or even larger than the pN contribution. 
These preliminary considerations showed that there might exist close compact binary 
systems for which the finite size effects introduced by rotational deformation are in 
the same range as the first post-Newtonian correction. Employing this assumption 
we shall now derive the quasi-Keplerian parametrization at the first post-Newtonian 
order including leading order quadrupole coupling (further on denoted as g-coupling). 
Introducing the reduced energy and angular momentum, E =: E/jji and J =: J / 
respectively, the IpN conserved energy including Newtonian g-coupling reads 



2 



GM 



la 



GM 

2r 



(3 + v)v'^ + vf' 



GM \ 
r J 



(7) 



Since the spin is parallel to the orbital angular momentum J the orbital plane is invariant 
in space. We are thus encouraged to introduce polar coordinates r and ip in the usual 
way. Inserting 
exact EoM to be 



-|- r^0^ into Eq. ([7]) one finds, after a little algebra, the IpN 



r 



E 



GM 

c^r 



(4 - 2u) 



. 2fi G 



D 



(8) 
(9) 



where 



A=2E 



1 + |(3.-1)| 

l + (7z/-6)- 



B = GM 

G = -J^ + ^ [2(1 - 3u)EJ^ 

G KA 1^ 

D = GMq\{%-'ivp^f- 



(5z/ - 10)G2A^'] 



(10) 



are constants. In the standard approach of Damour and Deruelle it is crucial that D 
is of order 0(c^^) and thus a small quantity. Now in our case D depends not only on 
c but it is also linear in the deformation parameter q. If the spinning component is 
governed by a soft EoS the correction to the orbital motion induced by the g-coupling 
is much larger than the IpN corrections. This is usual the case for main-sequence star 
binaries (see e.g. Claret and Willems 0). For compact stars in close binary systems, 
on the other hand, the contribution of the g-coupling can be of the same order as the 
IpN orbital correction. Under this assumption we can apply Damour and Deruelles 
strategy straightforwardly, deriving a quasi-Keplerian solution up to linear order of q. 
This yields 

sinw = n(i(; — to), (11) 



Cr cos u 



u 



if = 2(k + l)arctan 



'1 + e, 
1 - e,, 



■ tan ■ 



u 



(12) 
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where and depend on the coefficients defined above as 



B 
A 



1 



C 

2J2' 
DA 



GM D 



and n — ^ —A? jB"^, and k is given by 



et 



1 - 



A 
B^ 



C + 



BD 



K — 



J2 



+ 



2J2 



(13) 
(14) 

(15) 



As it turns out, the parameter which has to be small for this solution to hold is 
5 = g/j2(xl/c2. Expressing the parameters of the quasi-Keplerian solution in terms of 
the IpN conserved energy and 5, we find 



Or 



GM 

' 2E 

8E^ 



1 CM, 



1 - 



et = er 



1 + EU + 



-3iy 



2.1. Hamiltonian formulation 

The quasi-Keplerian solution given above describes the dynamics of a spinning compact 
binary with Newtonian quadrupole coupling at the first post-Newtonian order. However, 
according to GR the system looses energy due to the emission of gravitational waves, 
beginning at the order {v/cY in post-Newtonian approximation schemes. There are 
basically two ways to study the dynamical evolution of the binary system. In one 
approach the gravitational wave emission is considered as a secular effect and the 
dissipative terms do not enter the EoM. This allows the derivation of the quasi- 
Keplerian parametrization given above. Now we shall follow the other approach, where 
the radiation reaction terms are included into the EoM. In a first step we derive the 
Hamiltonian formulation to the conservative system. Energy loss due to emission of 
gravitational waves is incorporated by a time dependent radiation reaction Hamiltonian 
in a second step. 

The total Hamiltonian of a spinning compact binary system up to first post-Newtonian 
approximation reads 

HlpN = Hn + Hq + HpN, (16) 

where Hn and Hpjsf denote the point-particle Hamiltonian at Newtonian and first post- 
Newtonian order, respectively. If we would extend our analysis up to 2pN, other spin- 
depended terms would be present: at 1.5pN order the relativistic spin-orbit coupling 
enters into the Hamiltonian. It is, among others, responsible for the Lense-Thirring 
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effect. At the level of the second post-Newtonian approximation the relativistic spin- 
spin coupling leads to a precession of the orbital plane. Both, incorporating spin-orbit 
as well as investigating the spin-spin coupling, is beyond this letter, which is devoted to 
the study of the influence of certain finite size effects on the dynamical evolution and 
gravitational wave emission of the binary system. 

In the center of mass system the IpN Hamiltonian, including g-coupling, reads 



H 



IpN 



1 

2^ 



_GM 
2r~ 



GMfi 



r 



GMfi 1 

2r3 



3z/ - 1 



2 4 



2 2 

(3 + 2z/)^ + (3 + i/)^ 



+ 



2^2 



;i7) 



Note, that Hip^ = £ = fiE, since Hip^ is conserved at the first post-Newtonian order. 
As has been already mentioned before, the quadrupole interaction term is de facto a 
Newtonian correction to the point particle Hamiltonian. This is important to keep 
in mind when, for instance, calculating the orbital evolution and gravitational wave 
emission of binary pulsars with a main sequence star companion, such as PSR B1259-63 
[12] . For the Hamiltonian equations that govern the time evolution of the binary system 
one finds 



Pr 



Pr 

Pv 
fir"^ 

pI 



1 + 



3z/ 



1 + 



2 
3z/ 



Pr + 



PZ 



(3 + 2v) 



GMfi 



fir-^ r'^ 
+ ^ 



1 + ^ 
2r2 



GMii^ 



(3 + z/) 



+ - 



?>y-l Pt, 

2fi^ r^ 



2^p; 



GM 
2/ir2 



Pt 



(3 + 2z/K + 3(3 + z/)-f 



0. 



Now let us include leading order dissipative effects into that scheme. In general, the 
leading order energy dissipation of a matter distribution is governed by the time- 
dependent radiation reaction Hamiltonian [18] 



Hreac 



2G, 
5c5 



+ 



(19) 



p AtxG 

where is the symmetric tracefree mass quadrupole tensor of the matter distribution, p 
is the coordinate rest-mass density, vTj the momentum density and U is the gravitational 
potential that satisfies the Poisson equation with source term p. Since we are treating the 
g-dependent terms as being formally of first post-Newtonian order, the components of 
the binary can be considered as pointlike objects troughout the calculation of Hreacit). 
Thus in the center of mass frame the radiation reaction Hamiltonian describing the 
leading order energy dissipation due to gravitational wave emission is given by 

it) PlEl-GMp^ . (20) 



Hreacit^ 



2G (3) 
5c5 
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It is crucial to consider Xlj [t) as a function of time, and not as a function of generalized 
coordinates and momenta, when calculating the radiation reaction part of the EoM 
according to 

/ . \ dHfeac f 'i\ dHre.ar. 

yPiJreac r% • ; \Q )reac 



dpi 



(3) 

Only afterwards I>j can be expressed as a function oi pr,p^,r and (p. Explicitly, the 
calculation yields [1],[7] 

iPr^rad 
i,Pip)rad 
{r)rad 



8GV (GM'^v p 
3 r^c^ 



(0) 



rad 



\ 

8G^ (2GM 

5 vr^d' 
8 

15 z/r^c^ 
8 G^prP^ 



<fi 



3, ,2 



- 



pI 



2pI + 6% 



3 i^r^c^ 

For numerical calculations it is useful to introduce scaled variables such that G = c = 
The corresponding scaling is given by 

G^M^ 



(21) 
= 1. 



Pr = fiCPr 



Pv 



GMfi 

P^, 



GM 
— T^r, 



H = fic H, q 



Applying this, the Hamiltonian equations governing the evolution of the binary system 
including leading order radiation back reaction read 



Pr 



1 + 



3z/ 



P^ 
Pi 



1 + 



2 
3u 



3 + 2u 



r 

3 + u 



15 



2p? + 6 J 



r 



WPrP^ 



(22) 
(23) 



Pr ~9 



1 + ^' 
2f2 



3p-IpI 



5 



1 8z/ »r 

+ — H 

^3 3 



1 + 4 



f 



3 + 3 



(24) 



(25) 



Neither the total energy nor the orbital angular momentum is conserved, as indicated 
byEq. ([25]). 

The time evolution of binary systems described by Eqs. fl22|) -f l25|) is fully determined by 
3 parameters: the semi-major axis ar, the orbital eccentricity and the deformation 
parameter q, which have to be known for t = 0. Starting the numerical integration in 
the periastron, i.e. at ^^(O) = 0, the initial values for r and pr follow immediately as 

r(0)=ro = a,(0)(l-e,(0)), Pr(0) = 0. 

To determine the initial value for we use that, at the beginning of the integration, 
the total energy of the system is given by the conservative part of the Hamiltonian, or. 
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using the reduced energy, E{Q) = E = Hip]\i/fi. It follows then from Eq. (fT7|l thal| 



[{I - 3u)rl{E 



orb\2 
N ) 



+4(1 - iy)GMroE"^' + (6 - iy)G''M''] 



where E'^^ is the Newtonian energy of the orbit. 



3. Higher order gravitational wave emission 

In the previous sections we have investigated the dynamical evolution of a spinning 
compact binary system at the first post Newtonian approximation. Now we shall turn 
our attention to the gravitational waves emitted by the system. Far away from the 
source the space time can be assumed to be asymptotically flat, such that the metric is 
locally Minkowskian. In fact, in asymptotically fiat space-times the gravitational waves 
emitted by an isolated binary systems are expected to obey a multipole expansion of 
the form (see e.g. [T9] ) 

1-2 



TT 



G 



+ 



P, 



ijkm 



WE 



i-i 



1=2 

81 



(/ + 1)! 



4 
U 

m 



j^') (t 



D/c)Na,_ 



{t - D/c)N,Na,_ 



(26) 



the 



where h^^" is the symmetric-tracefree (STF) part of the metric perturbation /ijj, 
brackets denote symmetrization and D is the source-observer distance. The unit vector 
N points from the binary to the observer and Ai = aia2 . . . ai {ai = 1,2,3) is a multi- 
index, and J'ai are the STF mass and current multipole moments that parameterize 
the radiation field in a Cartesian coordinate frame. However, if the direction of the 
angular momentum is conserved, it is more suitable to use STF-multipole moments 



and S' 



Im 



m 



-I 



, /, that are irreducibly defined with respect to the axis of angular 



momentum. The relation of these to the Cartesian multipole components are given 
in Eqs. (1A.1|) and (]A.2|) in the appendix. From Eq. (!26|) it is then derived that the 
radiation field hj^ , expressed in terms of time derivatives of J'™" and S*'™, is given by 

G 



h 



TT 



EE 



1=2 m= 



-I 



+ I ^ 



1 

C 
l-l 



1-2 



lW"^{t-D/c)TP''"\Q,^) 



D/c)7;f''"(e,<i>) 



(27) 



where T^^^'^™ and T-^^'^™ are the so-called pure-spin tensor-spherical harmonics of electric 
and magnetic type. These harmonics are orthonormal on the unit sphere. In fact, 
introducing unit vectors $ and G, they can be decomposed into a term proportional to 
(6 (g) 6 - <l (g) $) and (0 (g) $ + $ (g) 0), respectively. That way, if the T^/^^M ^re 



% Note that fip^ and pr 



^ipr 



Next to leading order gravitational wave emission 



11 



known, one obtains the polarization states /i+ and hy^ oi the radiation field from Eq. 



fl27j) without any further calculations. The pure-spin tensor-spherical harmonics needed 
here are given by Eqs. ( lA.7[) -( |ArT6l) in the appendij|3- 

In section 2 the point particle contribution was taken into account up to first post- 
Newtonian approximation, and the quadrupole coupling term, though present already 
at Newtonian order, was assumed to be of the same order as the IpN corrections 
to the point particle dynamics. That means, we have to go beyond the leading 
order gravitational wave formula. Considering the dynamics up to IpN requires the 
application of the multipole expansion (!27|) up to / = 4 for the mass multipole moments 
and up to Z = 3 for the current multipole moments. Explicitly, neglecting all higher 
order terms, Eq. fl27|) is reduced tcEl 



G 



2 [-3 2 

j{2)2mrpE2,2m -\- - ) j{3)3mrpE2,3Tn _|_ g{2)2mrp 

m=— 2 1."^=— 3 m=—2 

r 4 3 

J j{4:)4:mrpE2,Am _|_ g{3)3mrp 



B2,2m 
X 



nB2,3m 
,x 



or, in a more convenient form. 



(1) 



,(2) 



+,x 



(2^ 



(29) 



c 

Note that according to Eq. (l28ll except for /(2)2m q]\ other time derivatives and multipole 
moments are required only at leading order. 

Let us start by noting that, in the center-of-mass system, the mass and current multipole 
moments of a two-body system are given by 

•2 _ 5-8uGM ' 

^2 



T. 



+ 



//(I — 3z/ 



1 + 3.)"- 

42 ^ 



7c2 



I, 



ijk 



I, 



ijkl 



2W 

-/iA/1 - Ai'X{iXjXk) , 

- ?,l')X(^iXjXkXl), 



12(r ■ v)xi^iVj) + llr^V{iVj) 



(30) 

(31) 
(32) 
(33) 
(34) 



Jij = -/iVl - 4:Ueab{iXj)XaVb, 
Jijk = yU(l - ^^)Bab(kXiXj)XaVb, 

where brackets denote the STF-part of the correspondent tensor (see appendix A). 
At this point it is neccessary to consider a moment the contribution of the stellar 
mass-quadrupole moment. Since one component of the binary is spinning and thus 
automatically gains a finite size there could be, in principle, a contribution of X?j to the 
gravitational wave emission of the system. This contribution is, however, very small, 
unless the energy stored in the internal stellar degrees of freedom, e.g. oscillations of 
the star, is comparable to the orbital energy. From now on we shall assume that I^j is 

+ In this representation, /i+ is the (8(g)8-$(Ki$)-part of Eq. ([77)1 . while fey is the (9(8)<i> + <l>(X)8)-part. 
* From now on we omit the supscript TT for notational convenience. 
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either trivial or can be neglected compared to all other terms present^ in the calculation. 
Using the Newtonian equation of motion for a point-particle binarj 
GM 



V 



to calculate j(3)3m^ j(4)4m^ jr(2)2m j{3)3fn ^^^j ^j^g equatiou of motion in the 
form 



GM 



V 



r 1 + 



3q 



+ 



GM 



2r2 

+ (4 - 2u){r ■ v)v 

-(2) 



r GM 



V 



Sv (r ■ v) 



to calculate X^- . one finds, after some lengthy calculations (see also [20] ) 



9 v'^ 



54z/ - 25 GM 



21 



25 + 9z/ GM , , GM 
+2iiXi,Vj) -^(r ■ v) - f^x^iXj)^^ 



2 , Ar-v 



2 



2 1 
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while the time derivatives of the current multipole moments read 
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In particular, within this approximation the q coupling is only relevant for the second 
time derivative of the mass quadrupole tensor. Using the relation between the two classes 
of multipole moments given by Eqs. (lA.ip and (lA.2p we get the following expressions 
for the time derivatives of P"^ and 5*'™, now in polar coordinates: 
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Remember that q is treated formally as a IpN quantity. 
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Here we used that VI - = (mi - m2)/M. Since only 5'^^)^^ I^^^^^ and /^^^^^ 
contribute to the first correction to the leading order quadrupole formula h^+]xi 
becomes clear that h^^]x vanishes for equal mass binary systems. In that case the first 
nontrivial correction to /i^^^^ is of order 1/c^. 

It is possible to derive, after some straightforward but rather lengthy calculations, 
analytic expressions for the time derivatives of the multipole moments in terms of e,.. E 
and 5; the corresponding relations are given in appendix C. For our purposes it is, 
however, more useful to express the polarization states of the gravitational radiation 
field in terms of generalized coordinates and velocities or - in Hamiltonian formulation 
- in terms of generalized coordinates and momenta. Inserting above relations in Eq. 
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Note that, in order to emphasize the character of the g-couphng, the g-dependent terms 
have been included into the leading order component of the radiation field. Defining 
Am = mi — 1712 the first correction terms read 
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Note that these expressions depend on the mass difference and vanish for equal mass 
binaries. The next corrections to the gravitational waveforms read 
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The polarization states of the gravitational radiation field, expressed in terms of 
generalized coordinates and momenta, can be found in appendix B. 



(5 



4. Discussion 



It has been long known that finize size effects introduce a periastron shift already at the 
level of Newtonian theory. For a couple of main sequence star binaries the total apsidal 
motion (ptot has been determined from observational evidence. Compared with the 
contribution iprei predicted by GR it became obvious that in all systems the Newtonian 
perturbations give the dominant contribution to iptot (for an overview see e.g. [9]). This 
is due to the " soft" equations of state governing the stellar matter of main sequence stars. 
For compact star binaries Newtonian perturbations are often neglected. In particular, 
it is often argued that the effect of the spin-induced quadrupole is too small unless 
the compact star (e.g. a neutron star) is rotating near the mass-shedding limit [3]. 
However, even for NS-NS or NS-BH binaries this argument does not hold completely. It 
has been shown in previous sections that for close NS-NS binaries the potential energy 
introduced by the coupling can be considerably larger than the corresponding 1.5pN 
orbital correction terms, though it is smaller by a factor 100 or more than the IpN 
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contribution. Thus, already at the level of the first post-Newtonian approximation in 
the EoM the rotational deformation induces a non-relativistic periastron shift, which 
accumulates over a large number of periods (figures [3]) andH]). In order to obtain highly 
accurate templates it is thus desirable to take into account these corrections properly, 
at least for fast spinning neutron stars in a compact binary system. 
With the space-bound laser interferometric detector LISA at hand the frequency band 
accessible to observations will be extended to much lower frequencies (10^^ to 10~^Hz), 
which enlarges the number of possible sources enormeously. In particular, with LISA 
not only BH-BH, NS-NS and NS-BH binaries should be detectable, but also white dwarf 
binaries. In particular to this class of compact binaries the analysis shown in this paper 
applies. It has been argued by Willems et al. in a recent paper |21] that - contrary to 
prevailing opinions - there might exist a class of eccentric galactic double white dwarfs, 
which are formed by interactions in tidal clusters. Willems et al. showed that tides 
and stellar rotation strongly dominate the periastron shift at orbital frequencies >1 
mHz. The phase shifts induced by these Newtonian perturbations are much larger than 
the general relativistic corrections then. They conclude that it is essential to include 
phase shifts generate by Newtonian perturbations into the signal templates in order to 
not bias LISA surveys against eccentric double white dwarfs. Generally, neglecting the 
contribution to (ptot induced by rotational deformation will lead to an overestimation of 
the total mass derived from (ptot- 

In this paper the competing influences of the rotational deformation and the IpN 
correction terms were examined in more detail. In particular, we succeeded in calculating 
a IpN quasi-Keplerian solution, which takes into account finite size effects up to linear 
order in the quadrupole deformation parameter q. The results given in section 2 are valid 
as long as q/ J"^ is of the order 0(c~^). For white dwarf binaries or binary pulsars such 
as PSR 1259-63 the periastron shift induced by rotational deformation is possibly much 
larger than the general relativistic contribution. In that case, Eqs. ( ITTlfT2l) still apply 
in the limit f/c — 0. In section 3 the polarization states of the gravitational radiation 
field are calculated beyond the leading order approximation. For non-spinning compact 
binaries the corresponding waveforms are shown in figures [1] and [2l In these figures 
waveforms calculated using the leading order expressions x'-"^ are compared to the 
IpN correct waveforms with the next to leading order corrections h^^^ and h^^y^ taken 
into account. The first correction, /i+''x, is nontrivial only for different mass binaries, i.e. 
for equal-mass binaries the first non-vanishing correction to the leading order formula 
appears at the order 0(c~^). 

The influence of the g-coupling on the gravitational waveforms is shown in figures [3] and 
m As expected, the spin-induced quadrupole moment leads to a phase shift compared 
to the pure point-particle GW emission. Moreover, the quadrupole deformation of the 
spinning compact objects speeds up the inspiral process, as has been shown in figures [5] 
and [6] for an equal-mass binary in a slightly elliptic orbit. 

More analysis is needed in order to fully understand the imprint of finite size effects onto 
the gravitational wave pattern of close compact binary systems beyond the leading order. 
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In particular it would be highly desirable to include the stellar oscillation modes into 
the calculations. From previous works it is expected that in these cases so called tidal 
resonances will have an important impact on the inspiral process and the gravitational 
wave emission of the binary [1], [5], |7j|. 



Acknowledgements 

I am grateful to Gerhard Schafer for helpful discussions and careful reading of the 
manuscript. This work is supported by the Deutsche Forschungsgemeinschaft (DFG) 
through SFB /TR7 " Gravitationswellenastronomie" . 



Appendix A. Useful relations 

The mass and current multipole moments J'™ and S*'™ (m = — /, ...,/) that are 
irreducibly defined with respect to the orbital angular momentum axis are related to 
and J7a, according to 



ji^,.. _ IStt / (/ + !)(/ + 2) 



^ ^^^""(/ + 1)(2/ + 1)!!V 2(/-l)/ ' ^^-^^ 



where, for m > 
and 



t = (-ir (2/ - 1)11 J + 4,) . . . + ,^?JC„ . . .^?„<A.3) 



rj™ = for m < 0. (A.4) 

The complex conjugates are given by 

Jim* = (-1)"^/'-"^, S^"^* = (-1)™^'-™. (A.5) 

The pure-spin tensor-spherical harmonics are orthonormal on the unit sphere. For the 
complex conjugate the following relation holds: 

rpE/B2,lm* ^ -^^mrpE / B2,l—m 

Defining 

T+ = e(g)0-<l>(g)<l>, T_ = e(g)<l + <l(8)0 

the expressions needed in the paper read 
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(A.9) 
(A.IO) 
(A.ll) 
(A.12) 
(A.13) 
(A. 14) 
(A.15) 
(A.16) 



Appendix A. 1. Symmetric tracefree tensors 



Throughout this paper symmetric-tracefree 3rd and Ath rank tensors are used. 
Symmetrizing a tensor of rank p requires to take the properly weighted sum over all 
index permutations, 

permut. 

The tracefree part of the tensor Tj^,, is calculated according to [22] 
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In particular, 3r(i and Ath rank STF tensors are given by 

T{abc) = T(^abc) — [SabT(cii) + hcT[aii) + ^acT{bii)\ , 
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Appendix B. Expressions for /i+ and vn terms of generalized coordinates 
and momenta 

The expressions for the leading and next to leading order contribution to the polarization 
states of the gravitational wave field, ^+ and h^, read 
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Appendix C. Higher order gravitational wave forms: Analytic expressions 



Using the quasi-Keplerian parametrization derived in section 2 it is possible to calculate 
analytic expressions for the times derivatives of the STF multipole moments entering in 
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the multipole expansion (1271) . Defining F{u) = 1 — cosm one obtains 
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Figure 2. /ix -component of a non-spinning, equal mass binary system. All parameters 
are the same as in Fig. [TJ 
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Figure 3. Influence of the quadrupole coupling on the gravitational wave emission. 
The IpN correct /i+-component emitted by a non-spinning binary with semi-major axis 
and eccentricity = 40 and = 0.3, respectively, is compared with the corresponding 
waveform emitted by a NS-NS binary with g = 4. The masses are mi = 3m2. Observer- 
dependent angles are $ = 7r/2,6 = 7r/4. 
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Figure 4. Influence of the quadrupole coupling on the gravitational wave emission. 
The IpN correct /i+-componcnt emitted by a non-spinning, equal mass binary with 
semi-major axis and eccentricity cir — 40 and — 0.6, respectively, is compared 
with the corresponding waveform emitted by a NS-NS binary with g = 4. Observer- 
dependent angles are $ = 7r/2, = 7r/4. 
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Figure 5. /i.-|--component of the gravitational wave field emitted by a non-spinning, 
equal mass binary during the inspiral process (initial values ar(0) = 50,6^(0) = 0.3). 
Observer-dependent parameters $ = 7r/2,6 = 7r/4. 
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Figure 6. /i+-component of the gravitational wave field emitted by an equal mass 
binary with q = A during the inspiral process (initial values 0^(0) = 50, 6^(0) = 0.3). 
Observer-dependent parameters $ = 7r/2, 6 = 7r/4. 



